Rockefeller University Report: RU 97/E-26 



o 



o 



Comments on the Erhan-Schlein model of damping the 
pomeron flux at small x-pomeron 

K. GOULIANOS 
(dino@pliysics.rockefeller.edu) 

The Rockefeller University, 1230 York Avenue 
New York, NY 10021, USA 

February 7, 2008 
ABSTRACT 

We explore the theoretical and experimental consequences of a model proposed by 
Samim Erhan and Peter Schlein for unitarizing the diffractive amplitude by damping 
the pomeron flux at small x-pomeron and conclude that the model is unphysical and 
contradicts well established experimental data. 

1 Introduction 

In Regge theory, the high energy behaviour of hadronic cross sections is dominated by 
pomeron exchange 0, |^. For a simple pomeron pole, the pp elastic, total and single 
diffraction dissociation cross sections are given by 
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where a{t) = 1 + e + a't is the pomeron Regge trajectory, Pippit) is the coupling of the 
pomeron to the proton, g{t) the triple-pomeron coupling, s' the IP — p center of mass 
energy squared, ^ = Xp = s'/s = M'^/s the fraction of the momentum of the proton 
carried by the pomeron, M the diffractive mass and so an energy scale not specified by 
the theory. 

In analogy with Eq. (Q), the term in brackets in (|^) is identified as the IP — p total 
cross section, and therefore the factor 
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is interpreted as a "pomeron flux" (per proton) and used in calculating hard processes 
in diffraction dissociation in the Ingelman-Schlein model 0. 

Experimentally, the triple-pomeron coupling was found not to depend on t [|T| and 
therefore we will use g{t) = g{0) and a^^{s',t) = aJ^^{s'). 

The function F{t) represents the proton form factor. Donnachie and Landshoff pro- 
posed [§| that the appropriate form factor for pp elastic and diffractive scattering is the 
isoscalar form factor measured in electron-nucleon scattering, namely 
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where m is the mass of the proton. When using this form factor, the pomeron flux is 
referred to as the Donnachie-Landshoff (DL) flux.Q 

In terms of the diffractive mass, the diffraction dissociation cross section is given by 
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where aj^^ = Prp^O) g{0)/ Sq. The total diffractive cross section is given by 

^ rO.ls roo 1 / s \ 2"'* 



where the lower limit of the integration is the effective diffractive threshold, Mq =1.5 
GeV^ and the upper limit, M^„^ = 0.1s, corresponding to ^ = 0.1, is dictated by 
the coherence condition for diffraction [|l[ (the contribution of the > 0.1s region to 
the integrated cross section would be, in any case, negligibly small). 

The ~ s^" dependence of asd{s) eventually leads to a diffractive cross section larger 
than the total, and therefore to violation of unitarity. The unitarity problem is a more 
general problem in pomeron pole dominance. It is well known that the power law s- 
dependence of the total cross section, ~ s^, violates the unitarity based Froissart bound, 
which predicts that the total cross section cannot rise faster than ~ In^ s. Unitarity is 
also violated by the s-dependence of the ratio aei/(TT ~ s% which eventually exceeds the 
black disc bound of one half {aei < \o't)- However, in both the elastic and total cross 
sections unitarization can be achieved by eikonalizing the elastic amplitude H]. 



2 The renormalized pomeron flux model 

Attempts to unitarize the diffractive amplitude by eikonalization or by including 
cuts have met with only moderate success. Through such efforts it became clear that 
these "screening corrections" affect mainly the normalization of the diffractive amplitude, 
but leave the dependence almost unchanged. Such a trend is clearly present in the 
data, as demonstrated by comparing the CDF diffractive differential pp cross sections 

^ The factor K in the DL flux is Kdl = {3f3]pq)'^ /An'^ , where is the pomeron-quark coupUng. 
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at y/s =546 and 1800 GeV with pp cross sections at y/s = 20 GeV. Motivated by these 
theoretical results and the trend observed in the data, a phenomenological approach 
to unitarizing the diffractive amplitude was proposed based on "renormalizing" the 
pomeron flux according to 

fNi^,t) = D-fjp/^i^,t) (8) 
where the factor D is determined by setting 
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D / fjp/pi^,mdt = i (9) 

if the value of the integrated standard flux exceeds unity (the limits ^min and 0.1 of 
the ^-integration are related to the limits Mq and 0.1s of the M^-integration in Eq. |^ 
by ^ = M^/s). Such a normalization, which corresponds to at most one pomeron per 
proton, leads to interpreting the pomeron flux as a probability density simply describing 
the ^ and t distributions of the exchanged pomeron in a diffractive process. 
The integrated standard flux is given by 
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and therefore D ~ s~'^^. Thus, flux renormalization approximately cancels the s^*^ de- 
pendence in Eq. (^) resulting in a slower rise of the diffractive relative to the total 
cross section with energy and preserving unitarity. Figure 1, which is an updated ver- 
sion (including more data) of a figure presented in Ref. [§, shows the total diffractive 
cross section as a function of y/s along with the predictions obtained using the standard 
(dashed curve) and renormalized (solid curve) pomeron flux (the dashed-dotted curve is 
the prediction of the Erhan-Schlein model, which is discussed below). Renormalized flux 
predictions of differential cross sections also show good agreement with data |jlO|, pT| . 
Finally, using the Ingelman-Schlein model the renormalized pomeron flux predicts 
correctly the measured rates of hard diffractive processes [|, 1^, . 

Predictions of hard diffraction rates using the standard/DL flux are unreliable due 
to a theoretical uncertainty inherent in the flux normalization. In Eq. (^ it is seen that 
/3jpp(0) can only be determined from the experimentally measured total cross section in 
terms of the energy scale Sq, which, as mentioned above, is not given by the theory (sq 
is usually set to 1 GeV^, the hadron mass scale, but this is only a convention). Thus, 
the normalization of the standard flux is unknown and therefore only predictions for 
relative hard diffraction rates are possible, as for example for the ratio of diffractive dijet 
production at two different energies. 

The normalization uncertainty is resolved in the flux renormalization model. The en- 
ergy scale. So, can be determined by setting the flux integral to unity at ^/s ~ 20 GeV, 
where the total diffractive cross section turns from rising as ~ s^^ to assuming a rather 
flat s-dependence, presumably due to the saturation of the pomeron flux (see Fig. 1). 
For ^/s > 20 GeV, where the flux integral is unity, the normalization is self-determined 
from Eq. (P). Thus, predictions for hard diffraction using the renormalized pomeron 
flux can be made not only for relative but also for absolute rates. The recently reported 
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CDF diffractive W and dijet production rates [jl2|, [T^ are in excellent agreement with 
the renormalized flux predictions |10|, |TT|] . 



3 The Erhan-Schlein model 



Erhan and Schlein have taken a different approach to solving the unitarity problem 
of the triple-pomeron amplitude |1^]. They introduce a factor D{C,) in the flux to 



damp the small-^ values and thus slow down the ~ s dependence of the flux integral. 
The normalization of the flux is left unchanged. To slow down the rise with s of the 
differential cross section at the higher ^-values that are not affected by the damping 
factor, they introduce a IPIPR term, whose contribution increases at low energies due 
to the ~ dependence of the i?-term. To fit ISR and UA8 data at |)f:| ~ 1 — 2 GeV^, 

two more parameters are introduced. The detailed form of the proposed diffractive 
differential cross section is the following: 
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Using Eqs. ([TT]-|T5D and the fixed parameters 

e = 0.115 A = 0.0^ 



(11) 
(12) 

(13) 
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diffractive cross section data are fitted with the following four free parameters determined 
from the data: 



b 
a" 



0.73 ± 0.09 mb GeV-^ 
0.75 ± 0.27 GeV-2 
0.075 ± 0.017 GeV-^ 
5.0 ±0.6 



(17) 
(18) 
(19) 
(20) 



It is claimed that good agreement is obtained with "all available data" |]TB[. Below, we 
comment on the effect of the extra parameters introduced in the standard triple-pomeron 
amplitude, compare the predictions of the model with experimental pp and pp data and 
highlight some of its experimental and theoretical implications. 
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3.1 The parameters of the Erhan-Schlein model 

In addition to the introduction of the IPIPR term and the parameters associated with 
it, several additional parameters are introduced to the standard triple-regge amplitude: 

• A different pomeron intercept is used in the (T^^(s') term than in the flux factor 
(see Eq. |16|). 

• Two damping factors, -D(^), are used for two different ^-ranges. Figure (2a) shows 
D(^) as a function of ^. The roller-coaster shape of this distribution is mapped 
into the double resonance like distribution of the production cross section for a 
fixed diffractive mass at t = as a function of energy, as shown in Fig. (2b). 

• The effect of the parameter h of the e^* term in Eq. (p!2| ) is shown by the solid 
line in Fig. (2c). In the region 1 < \t\ < 2 GeV^ of the UA8 data, which have 
been presented in comparisons with the predictions of this model [|1^, |I5[, the h- 
parameter accounts for a t-dependent reduction of the cross section by a factor of 
2-5. 

• The effect of the parameter a", which introduces a dependence of 1/.^^""*^, is shown 
in Figs. (2c, 2d). Figure (2c) displays this factor for ^ = 0.05, the average UA8 
^-value, as a function of t. Within the UA8 t-range, it varies from 1.5 to 6. The 
^-dependence as a function of t for this factor is shown for various t-values in 
Fig. (2d). Again, within the UA8 range of t and ^, the effect is seen to be quite 
large. 

3.2 Comparison with experimental pp and pp data 

In Fig. 1, which shows total pp and pp diffractive cross sections as a function of -y/i, the 
dashed line represents the standard flux prediction, the solid line the renormalized flux 
prediction and the dashed-dotted line the Erhan-Schlein prediction. The latter exibits 
two turn-overs as -^i increases, one at ~ 10 GeV and the other at ~ 200 GeV. The 
^min ~ 1/s values corresponding to these energies are at the interfaces of the damping 
factors that are used in the model. Clearly, a single damping factor would not provide 
a good fit to the data. 

Figure (3) shows differential cross sections ^(Pa/d^dt at t = —0.05 GeV^ as a function 
of ^ for fixed target pp data at ^/s = 14 and 20 GeV [jT^ and for the CDF pp data at 
546 and 1800 GeV p . The fixed target data are for masses above the resonance region 
and the CDF data are for ^-values large enough not to be affected by the experimental 
resolution of the ^-measurement. The dashed curves are fits to the data using the form 

with ajp(t) = 1.104 + 0.25t g] and a^{t) = 0.9t. The first term in Eq. (|T]) is the 
form of the triple pomeron amplitude and the second term the form for reggeized pion 
exchange. Fits with IPIPIP and vr-exchange terms have been shown to represent well the 
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pp/pp diffractive data [jl], |I0|, |ri|]. Using the renormalized flux, such fits with only one free 
parameter, namely the triple-pomeron coupling constant, yield differential cross sections 
that are in good agreement with the data not only in shape but also in normalization P, 



10, 11 



The solid curves in Fig. (3) were calculated using the Erhan-Schlein model. Three 
features of these curves are immediately apparent: 

• At low energies (Figs. 3a, 3b), the Erhan-Schlein cross sections are falling sharply 
with energy. This behaviour is due to the dependence of the (T^^(s') term in 

the IP PR amplitude (second term in Eq. |Tip, which dominates the low energy cross 
sections in this model due to the s~'^^'^^ factor, as compared to the dependence 
of the IP IP IP term (first term in Eq. |14|) . 



At the high energies (Figs. 3c, 3d), the curves initially rise as ^ decreases and then 
bend over and fall as ^ decreases further. The bending point is at ^ = 0.015, which 
is the ^-value where the -D(0 damping takes effect. 

At ^ = 0.035 the predictions are in relatively good agreement with the data. 
However, this agreement is not very meaningful as the parameters of the model 
were determined by fitting data of cross sections as a function of energy at this 



particular value of ^ [0, [T^. The argument offered for using only data at ^ = 0.035 
in the fits is that this ^-value is low enough for the data to be background-free, but 
also high enough so that the ^-distributions are not distorted by the experimental 
resolution in measuring ^. As mentioned above, however, none of the data points 
in Fig. (3) are affected by resolution and therefore the high energy data in the 
region of ^ < 0.035 provide a good testing ground for the model. Figures (3c, 3d) 
show that the model is not very successful in this region. 

3.3 Experimental and theoretical implications 

In addition to the experimental issues already discussed in connection with pp/pp cross 
sections, damping the pomeron flux at small-^ has serious implications for the HERA 
deep inelastic scattering (DIS) and photoproduction results. These results are in the 
range 10""^ < ^ < 10"^, within which the damping factor D{^) = 1 - 2700(^ - 0.015)^ is 



in effect. The DIS HI ||T^ and ZEUS [|18| results and the HI photoproduction results ||T9[ 
are compatible with a pomeron flux ~ 1/^^+^^ with e ~ 0.11. From ^ = 10~^ to ^ = 10^^ 
the damping factor D{C,) increases by a factor of 2.2 while with e from Eq. (plGj) 
decreases by a factor of 2.8. Therefore, if .^-damping were indeed in effect, the value of 
e expected at HERA would be e ~ 0.02. Thus, the HERA data contradict the small-^ 
damping hypothesis. 

On the theoretical side, there is no reason to expect that screening corrections should 
damp the cross section preferentially at small diffractive masses. In fact, as already 
mentioned, eikonalization leaves the distribution largely unchanged 0. As for the 
other parameters introduced in the model, there is no obvious reason why A should be 
different from e, or, if it were different, that it should have the value 0.08; neither is 
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there a reason why the term e should be needed in diffraction, since a corresponding 
term (e''*)^ = e^'^* is not needed in the form factor for elastic scattering. Finally, the 
introduction of the term a"t^ in the pomeron trajectory makes both the diffractive and 
the elastic cross sections blow up at large values of t. 



4 Conclusions 

In an effort to solve the unitarity problem inherent to the triple-pomeron description of 
the diffractive cross section, Erhan and Schlein introduce a ^-damping factor, D{^), that 
decreases the pomeron flux for ^ < 0.015. Since for ^ > 0.015 the flux is left unchanged, 
the differential cross section (Pa/d^dt given by the IP PIP amplitude still rises as ~ s'^ 
in this region, while experimentally it is found to decrease as ~ s~^. To balance the 
s-dependence, a IPIPR term that varies with s as ~ s~^'^^ is introduced in the model. 
However, as this term has a ^-dependence sharper than that of the IPIPIP term (a factor 
of is replaced by ^~°'^^), it is now more difficult to fit the ^-distributions of the data. 
To obtain better fits, more parameters are introduced. In addition to choosing A 7^ e, 
four free parameters are used. These do not include the three parameters needed for the 
damping factors (two ^-values and the value D = 0.4 at ^ = 10~^), which were chosen 
to make the model agree with the s-dependence of the total diffractive cross section. 
Including these three parameters and the parameter A in the list of free parameters, a 
total of 8 free parameters are used, not taking into account the shape of the damping 
factors, which was also chosen to optimize the fits. Furthermore, Erhan and Schlein 
point out that their model does not agree with data at >~ 0.5 GeV^, which includes 
the UA8 data, unless the ^-damping factor is not used in this region [|15|. However, 
no prescription is offered as to what happens at the interface between \t\ < 0.5, where 
damping is needed, and |t| > 0.5, where it is not. Presumably more free parameters, 
which are currently hidden, would be required to take this effect into account. Despite 
the large number of explicit and hidden free parameters, the predictions of the model 
are contradicted by both the pp/pp and the HERA data. The model is also theoretically 
unsound as it predicts the vanishing of the t = cross section for small and increasingly 
larger diffractive masses as the energy increases. 

In conclusion, the proposed model is ill-defined (low versus high \t\ behaviour), the- 
oretically inconsistent and unphysical, contradicts pp/pp and HERA experimental data, 
and has no predictive power for hard diffraction rates within the framework of the 
Ingelman-Schlein model. 
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Figure 1: The total pp/pp single diffraction dissociation cross section for ^ < 0.05 as a 
function of center of mass energy. The dashed curve is the triple-pomeron prediction, the 
sohd curve the renormahzed flux prediction and the dashed-dotted curve the prediction 
of the Erhan-Schlein model of damping the pomeron flux at small-^. The latter exibits 
two "turn-overs" as ^/s increases, one at ~ 10 GeV and the other at ~ 200 GeV. The 
Cmin ~ 1/s values corresponding to these energies are at the interfaces of the damping 
factors that are used in the model. 
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Figure 2: The parameters of the Erhan-Schlein model: (a) the damping factor, 
D{^), as a function of ^; (b) (Pa/dM^dt\t=o versus y/s for = 5 GeV^ (sohd) and 
10 X (fa/dM'^dt\t=o for = 50 GeV^ (dashed); (c) the factors e** (sohd) and T^""*' at 
^ — 0.05 (dashed) versus (d) the factor versus ^ for i = 0, 1, 1.5 and 2 GeV^. 
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Figure 3: Cross sections ^d?a /dtd^ versus ^ for pp/pp data compared to theoretical 
predictions. The dashed hues are fits with a triple-pomeron term and a pion exchange 
term (see Eq. ^T]). The sohd hues represent the predictions of the Erhan-Schlein modeL 
At the two highest energies, as ^ decreases the sohd curves are seen to "bend-over" and 
fall at ^ ~ 0.015, where the ^-damping factor becomes effective. 
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